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Abstract 

In this paper we consider a quantum open system and treat the master equation with 
some restricted dissipator which consists of a set of projection operators (projectors). The 
exact solution is given under the commutable approximation (in our terminology). This 
is the first step for constructing a general solution. 

In this paper we revisit dynamics of a quantum open system. First of all we explain our 
purpose in a short manner. See ^ as a general introduction to this subject. We consider 
a quantum open system S coupled to the environment E. Then the total system S + E is 
described by the Hamiltonian 

H s+E = H S ®1 E + 1 S ®H E + Hj 

where H 5, H E are respectively the Hamiltonians of the system and environment, and Hi is the 
Hamiltonian of the interaction. 
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Then under several assumptions (see [T]) the reduced dynamics of the system (which is not 
unitary !) is given by the Master Equation 

^P=-i[H a ,p]-V{p) (1) 
with the dissipator being the usual Lindblad form 

v ^ = \J2 { A ) A iP + p a ) a > - 2A iP A i) • ( 2 ) 

{i} 

Here p = p(t) is the density operator (matrix) of the system. 

It is not easy to solve the equation (DO) with the dissipator (T5]), so we make a simple and 
convenient assumption. Namely, the generators {A,} are given by Aj = ^f\jPj with projectors 
{Pj} ; pj = p h pj = p h p 3 p k = S jk P k . Note that we don't assume the rank Pj =1 (extended 
models). Then the dissipator becomes 

V (P) = \J2 X J ( P JP + P p i - 2P iP P i) ( 3 ) 

where {Aj} are decoherence parameters to determine the strength of the interaction. See [2], 
[3] (in [3] there is a very compact description on this subject). It is interesting to rewrite ([3]) as 

V (P) = \ll X j " P i) + P i)P P A = ll2 X i ( P iPQ + QiP P i) ■ ( 4 ) 

{j} {j} 

Note that {Qj} are also projectors satisfying PjQj = QjPj = for j G {j}. 
As a result we have only to solve the equation 

^ t P = ~< H P -P R )-\Y. X > ftPQi + ®iP P i) ( 5 ) 

{j} 

where we have set H = H$ for simplicity. 

In order to attack the equation (jSJ) let us make some mathematical preliminaries. For a 
matrix X = (x^) G M(n; C) we correspond to the vector X G C" 2 as 

X = (Xij) > X = (Xn, X12, • ■ • , X\ m , X n \, X n 2, • ■ ■ , x nn ) T (6) 

where T means the transpose. Then the following formula is well-known 

AXB = (A ® B T )X (7) 



for A,B,X£ M(n; C). Since the proof is easy we leave it to readers. 
By use of the formula the equation can be rewritten as 

j t p= |-i(ff®i-i®s T )-^A i (p i ®gj + g i ®p/ , )J^ ( 8 ) 

therefore the formal solution is given by 

pit) = exp | -it(H ®l-l®H T )-t ^(A,/2) (P j ® Qj + Q j ® Pf) | p(0). (9) 

To calculate exp(- ■ ■ ) explicitly is (almost) impossible, so we must appeal to some approxi- 
mation method. For that let us remind the Baker-Campbell-Hausdorff (B-C-H) formula. For 
A,Be M(n; C) we want to decompose as 

e A+B = eV^e 5 (10) 

The "interaction" term I(A, B) is given by 

I(A, B) = -l[A, B] + X - {[[A, B],B] + [A, [A,B]]}-\ . (11) 

The proof is easy. In fact, e^ A ' B ^ = e~ A e A+B e~ B by (jTUl) and we have only to apply the B-C-H 
formula ([I] and see also [5] as an interesting topic) 

e x e Y = e x+Y+(y2)[x,Y]+(y^){[[x,Y],Y]+[x,[x,Y]]}+- for X , Y G M(n; C) 

two times. 
For 

A = -it(H <g> 1 - 1 ® P T ), P = -t ^(A i /2) (P, ® Qj + ® P/) 

there is no method to calculate e I( - A ' B ^ explicitly as far as we know. Therefore we ignore this 
term, namely let us call it the "commutable approximation" . 

Under the commutable approximation we have only to calculate 
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pit) « exp {-it(H <8> 1 - 1 <g> P T )} exp i -t (P, g> + <g> Pf ) \ p(0) 



{./} 



(e^ ® e itHT ^j exp | -t \/2) {P j ® Qj + Qj <8> Pf ) j p(0) 

(e^) T ) exp | -t E( V 2 ) ® + ® *?) | P(0)- 



e ~HH o („itH\ T 



Next let us calculate the second term in (fl2l) . which is not so difficult as follows. 

(0) = exp J -t (p,- ® gj + g, ® pJ) 
{ {3} 

= H exp { (-a^/2) (p, ® gj + g, ® Pf) } 

{il 

= JJ { 1 ® 1 + (e-^*/ 2 - 1) (P, <g> gj + Qj <8> Pf) } 

0'} 

where we have used facts 

(a) {Pj ® gj + Qj ® Pf I j G {j}} are projectors commuting with each other. 

(b) e XR = 1 + (e A - l) R if P is a projector. 

Here we set Rj = Pj <8> gj + g^ ® Pf . For i < j < k we obtain 

(c) PiP, = (Pi ®Qf + Qi® pf ) ® gj + g, <8) pf ) = p ® pf + Pj ® pf . 

(d) R t RjR k = (P ® Pf + Pj ® Pf ) (P fc ® g^ + g fc (8) if) = 0. 
From ([IS]) and (c), (d) 

(tt)=n( i ® i +( e " v/2 - 1 )^} 

{i} 

=i®i+x; ( e_Ajt/2 - 1) ^ + e ( e_Ajt/2 - 1) ( e " Afc * /2 - 1) ^* 

j 3<k 

= 1 <8> 1 + ^ ( e ~ Xjt/2 - X ) ( P i + ^ P f ) + 

j 

£ (e-V/2 _ ^ ( e -A fc */2 _ 1) (p. g) pT + p fc ^ pT) 

i<fc 



Therefore 



p(t) *(e- itH ® (e^) T ) + £ (e~ A ^ 2 - l) (P, ® Qj + Q j ® if) + 

£ (e-A^/2 _ ^ ( e -A fet /2 _ ^ (p. g, pT + p fc s pT) I 



(15) 



(16) 



Coming back to matrix form by use of (JTj) we finally obtain 

p(t) « e"** jp(0) + £ (e^/ 2 - 1) (P i p(0)Q i + Q,p(0)P 

£ (e-M/a _ 1) (e-^'/ 2 - 1) (P,p(0)P fc + P fc p(0)P,) 1 e 
j<k ) 

or 

p(t) » e-** jp(0) + (e" Ajt/2 - 1) (PiPMQi + OiP(O)Pi) + 

I ^ ( e -M/a _ X ) ( e -A,t/2 _ x j (p. p(0 )p fe + P fcP (0)P,) 1 e UH (17) 
Z rtk J 

for j, k G {j}. This is the main result. 

A comment is in order. In the two qubit system a general density matrix is written as 

p(t) = - (l 2 ® 1 2 + Piitfdi <g> 1 2 + <?j(t)l 2 ® CTj + Tij(t)(Ti <S> (Tj) 

where we have used the Einstein's notation on summation . Using this expression one tries to 
solve the equation coming from pure decoherence term 

{j} {J} 
The equation is then reduced to a set of (relatively simple) equations of {p}, {q} and {r}. 

However, such a method (trial) is irrelevant as shown in the paper. Our method is quite 

general ! 

In this paper we considered the master equation with the dissipative being a set of projectors 
and constructed the exact solution under the commutable approximation. This is just the first 
step for constructing a general solution for the equation. 



In order to take one step forward we must take the "interaction" term I(A, B) in (TTTT) into 
consideration. However, such a method to calculate it has not been known as far as we know. 
Therefore it may be reasonable to restrict our target to some simple models. Further work will 
be needed and we will report it in a near future, [6] . 

On the other hand we are studying some related topics, see [7j and [8]. However, we make 
no comment on them in the paper. 

Lastly, we conclude the paper by stating our motivation. We are studying a quantum 
computation (computer) based on Cavity QED (see [S] and [ID]), so to construct a more re- 
alistic model of (robust) quantum computer we have to study a severe problem coming from 
decoherence. This is our future task. 
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